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ABSTRACT 

In this paper we discuss models of the X-rays and TeV 7-ray emission from 
BL Lac objects based on parallel electron-positron or electron-proton beams that 
form close to the central black hole owing to the strong electric fields generated by 
the accretion disk and possibly also by the black hole itself. Fitting the energy 
spectrum of the BL Lac object Mrk 501, we obtain tight constraints on the 
beam properties. Launching a sufficiently energetic beam requires rather strong 
magnetic fields close to the black hole (~ 100 — 1000 G). However, the model 
fits imply that the magnetic field in the emission region is only ~ 0.02 G. Thus, 
the particles are accelerated close to the black hole and propagate a considerable 
distance before instabilities trigger the dissipation of energy through synchrotron 
and self-Compton emission. We discuss various approaches to generate enough 
power to drive the jet and, at the same time, to accelerate particles to ~20 TeV 
energies. Although the parallel beam model has its own problems, it explains 
some of the long-standing problems that plague models based on Fermi-type 
particle acceleration, like the presence of a very high minimum Lorentz factor of 
accelerated particles. We conclude with a brief discussion of the implications of 
the model for the difference between the processes of jet formation in BL Lac 
type objects and in quasars. 

Subject headings: galaxies: jets — galaxies: BL Lacertae objects: individual 
(Mrk 501) — gamma rays: theory — X-rays: galaxies 
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Introduction 



1.1. Observations and models of the continuum emission from blazars 



Observations with the EGRET Energetic Gamma-Ray Experiment Telescope on board 
of the Compton Gamma-Ray Observatory (CGRO) revealed that blazars are powerful and 
variable emitters, not just at r adio through optical wavelengths but also at >100 MeV 7-ray 
energies (IHartman et al.lll999h . The 66 blazars detected with EGRET were mainly quasars, 
i.e. Flat Spectrum Radio Quasars and Optically Violent Variables. Observations with ground 
based Cerenkov telescopes showed that BL Lac objects, the low- power counterpart s to the 
high-powered EGRET quasars, emit even more energetic 7-rays (jPunch et al.lll992l ). Based 
on observations with ground based Cerenkov telescopes, more tha n a dozen BL Lac type 



objects have been iden tified as sources of >300 GeV gamma-rays (lAharonian et al.l 12005 
Horan fc Weekesll2004l ). 



The MeV and Te V 7-ray emission from blazars is commonly thought (e.g. iTavecchio 



2OO5I : iKrawczynskil l2006l ) to originate from r elativistic particle dominated outflows (jets) 



from mass accreting supermassive black holes flLynden-Belllll969l : IZeldovich &: Novikovlll971 



Reeslll978l ). The jet may form electromagnetically or through magnetohydrodyn amic pro 
Electromagne tic models come in two flavors. Either the a.ccreti on disk (ILovelace 



cesses. 



19761 : iBlandfordI Il976l ) or the Kerr black hole (IBlandford fc Znajekl 119771 ) launch a Poynt- 
ing flux dominated flow. The mechanism to convert the Poynting dominated outflow into 
a particle dominated one is not understood. In other models m agnetohydrodynamic pres 



sure may play a dominant role in the process of jet formation (IBlandford fc Paynd Il982 



Begelman. Blandford. fc Reed Il984j ). Models usually assume that the particle dominated 
outflows move with velocities v ~ c and Bulk Lorentz factors F = (1 — (f/c)^)~^/^ be- 
tween a few and ^50. Some models suggest very hig h bulk Lorentz factors between 100 and 



1000 flPohl fc Schhckeiserl I2OOOI : iRebillot et al.l 120061 ). At shocks within the jet plasma, the 
Fermi mechanism may transfer a fraction of the bul k kinetic en ergy into random kinetic en- 
ergy of high-energy particles, electrons, or protons (IReesI 119781 ). These accelerated particles 
themselves, or secondaries produced in cascades, may emit the observed continuum emission. 
Although the spectral energy distributions (SEDs) of blazars are often only sparsely sampled, 
most blazars seem to emit two distinct emission components, one at low energies and one 
at high energies. The two components are attributed either to the emission by a single par- 
ticle popul ation emitt i ng photons of vastly different energies through two different emission 
proce sses (iReed 119671 : IBlandford &: Reed 1 19781 : iKonigll Il98ll : iGhisellini. Maraschi. fc Treves 
19851 ). or to the emission by two particle populations. An example for the flrst are syn- 
chrotron Compton models in which a single population of electrons emit the low-energy and 
high-energy emission components as synchrotron and inverse Compton emission, respec- 
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tively. In synchrotron self-Compton (SSC) models, the synchrotron photons are the dom- 
inant source of target photons for inverse Compton processes. Examples for the latter are 
hadronic models in which the low-energy component originates as synchrotron radiation from 
a population of low-energy electrons. The high-e nergy component is synchrotron emission ei- 



ther from extreme ly high energy (EHE) protons (lAharonian Il2000l : iMiicke fc Protherodl2001 



Miicke et al.ll2003l ). or from seconda ry e"^ /e resulting from a synchrotron and pair-creation 
cascade initiated by EHE protons (IMannheinj Il993l) or high - energy electrons or photons 



(ILovelace. MacAuslan. fc Burns! Il979l : 



Levinson fc Blandford 



19951). 



Burns fc Lovelacdll982l : iBlandford fc Levinsonlll995 



Following the end of the CGRO in 2000, studies of the MeV emission from quasars 
has to await the launch of the next space-borne 7-ray telescope. The GLAST ( Gamma-Ray 
Large Space Telescope) satellite will be one order of magnitude more sensitive than EGRET 
and is scheduled for launch in the near future. Ground-based 7-ray observatories continue 
to provide data on the TeV 7-ray emission from BL Lac objects. Numerous broadband 
multiwavelength observation campaigns on a few objects (most notably Markarian (Mrk) 
421 (z = 0.031), Mrk 501 {z = 0.034), and lES 1959+650 {z = 0.048)) have yielded 
important results which shed some light on the emission mechanism. The key-results from 
these campaigns are: 



1. " 

( 


Tiere is a hig 


ily significant flux correlation between X-rays and TeV 7-rays 


Taka 


lashi et al. 


1996; 


Buckle V et al. 


1996; 


Krawczvnski. et al. 


2000; 


Sambruna et al. 



2000 



Fossati et al.l 120041 ). The "time lag" between the flux var iations at X-rays and 



at TeV 7-rays is s ufficiently short ( <1 hr) to evade detection (IMaraschi et al. 



1999 



Fossati et al.ll2004l ). The X-rays and TeV 7-rays are emitted close to the peaks of the 



low-energy and high-energy emission components, respectively. 

The sources show strong correlated X-ray/TeV 7-ray flares with X-ray and TeV 7-ray 
flux changes by factors of between a few and 20 on time scales between 15 min (Mrk 421) 
and a few hours (Mrk 501, lES 1959+650) (iGaidos et al. I Il996l : iKrawczynski. et al. 



2000). 



The campaigns did not reveal a highly signiflcant flux correlation of the radio to optical 
emission with the X -ray or TeV 7-ray emissions. The interested re ader is referred to 
Buckley et al. (Il996 ) f or wea k evidence for such a correlation and to Blazejowski et al. 



( 2OO5I ): Rebillot et al. (2006) for campaigns which did not reveal supportive evidence. 



4. 



As discussed further below, shock acceleration theories predict that the X-ray and 
TeV 7-ray energy spectra are softer during the early rising phases of flares than during 
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the later phases. While some flares showed such a behavio r, others did not (e.g 



Takahashi et al.ll2000l : iFalcone et al. II2004I : iFossati et al.ll2004f ). 



Despite some disparity of the data on a whole, the good X-ray/TeV-7-ray correlation 
strengthens the case for leptonic synchrotron-Compton models. 



1.2. The reference model and its problems 



In this paper, we focus the discussion on the X-ray and TeV-7-ray emission from BL 
Lac type objects. In Section HI we will briefly outline the relevance of the model for quasars. 

We use the term "reference model" for the following combination of model components. 
The accretion system launches a Poynting flux dominated jet and as the outflow propagates, 
the flow transforms into a particle dominated outflow. Shocks within the jet transfer a 
fraction of the jet's bulk kinetic energy to a few high-energy particles that emit synchrotron 
and Inverse Compton emission. 

The reference model suffers from a number of weak links. No simple mechanism has yet 
been suggested to convert the Poynting flux dominated outflow into a particle dominated 
outflow. Furthermore, the hypothesis of shock acceleration of electrons to ~TeV energies is 
observationally very poorly supported. The modeling of the SEDs of some BL Lacs require 
"non-standard" electron energy spectra. 

Models of Mrk 501 and Mrk 421 data require a minimum Lorentz fa ctor of ac 
celerated particles 7rnin on th e order of 10^ in the jet reference frame (e.g. iPian et al. 



19981 : iKrawczynski et al.l 120021). or non-power- law distributions with very high characteristic 
Loren tz factors fe.g. ISauge fc Henri 2004; Katarzyhski et al. 2006 : Giebels. Dubus fc Khelifi 
20061 ). If the shocks are internal to the jet and the jet medium is made of cold protons, simple 
arguments lead to 7min- values close to the proton-to-electron mass ratio 7min ~ trip /me = 
1836. If the shocks are external and the jet medium runs into a much slower target medium, 
the 7min-values could be higher by a factor of F. A modification of the reference model might 
be able to account for the high 7min -values and for the non-sta ndard electron spectral indices 
even in the internal shock model. iBykov fc Meszaros I (jl996[ ) point out that statistical ac- 
celeration by relativistic magnetohydrodynamic fluctuations in flow collision regions of jets 
might give rise to hard electron energy spectra and to 7min- values on the order of 10^. 

As mentioned above, the theory of shock acceleration predicts that the X-ray and TeV 
7-ray spectra of flares are soft during the early rising phases of flares. The beautiful Mrk 421 
observations taken in 2001 with the RXTE (Rossi X-ray Timing Explorer) show both, soft 
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and hard energy spectra during the early rising phases of flares ( iFossati et al.ll2004l . and Fos- 
sati and Buckley, private communication, 2006). This negative result can still be explained 
in the framework of the reference model, if a strong flare is made of the superposition of 
many small flares, or if other processes (e.g. the growth and decline of the shock, a changing 
viewing angle, radi ative cooling, and adiabati c particle losses) dominate the temporal evo- 
lution of the flares (Kirk &: Mastichiadisl ll999h. ISokolov. Marscher. fc McHardvl tooi ) and 
Sokolov fc Marscherl (120051 ) emphasize that the geometry and structure of the acceleration 
region also influence the observed spectral behavior. 

A different problem concerns the bulk Lorent z factors of the emitting jet plasmas. Sim- 



ple S S C models require bulk Lorentz factors F >25 (IKrawczynski et al.ll200ll : lKonopelko et al. 
20031 : iHenri &: Saugel 120061 ). While Very Large Baseline Array (VLBA) observations of 



quasars have rec ently succeeded in finding sources with app arent superluminal motions that 
support F ~ 50 ( Finer. Bhattarai. Edwards, fc Jonesll2006l ). the observations of BL Lac ob- 
jects in ge neral JListerl I2OO6 ) and the T eV sources Mrk 421, Mrk 501, lES 1959+650 in 



particular (IFiner fc Edwards! |2004| . |2005| ) show only subluminally moving components. In 
contrast to SSC models, "External Compton" models do not require such high bulk Lorentz 
factors. In these latter models an otherwise unobserved radiation component originating 
from a different emission zone than the X-ray and TeV 7-ray emission, provides the seed 
photons f or the inverse Compton processes that p rodu ce the observed 7-rays. Th e recent 



models of iGhisellini. Tavecchio. fc Chiabergd (120051 ) and lGeorganopoulos fc Kazanad (120031 ) 
can fit the data with lower bulk Lorentz factors. 



1.3. Structure of the paper 



Motivat ed by the probl e ms of the r e 



'erenc e mod el, we explore here models similar to 



the ones of iLovelacd (Il976l ). iBlandfordI (Il976l ). and iBlandford fc Znajekl (119771 ). As the 



models invoke parallel electron-positron or electron-proton beams, we will refer to them in 
the following as "parallel beam" models. A disk carrying a magnetic field and possibly 
the Kerr black hole induce a strong electric field parallel to the rotation axis of the black 
hole/disk system. We assume that the electric field will give rise to a beam of high-energy 
electrons and positrons or protons that move nearly parallel to the rotation axis of the 
accretion system. These particles, in turn, directly emit the observed X-ray and TeV 7-ray 
emission as synchrotron and inverse Compton emission, respectively. This model does not 
require the formation of shocks or Fermi-type acceleration mechanisms as in the reference 
model. In this paper, we scrutinize for the first time how the model can be applied to explain 
the X-ray/TeV 7-ray data that have recently become available. We will use the data to guide 
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the evaluation of the models. In Section [2l we will first infer the properties of the nearly 
parallel beam from modeling a specific data set. Once that we know the properties of the 
beam, we will discuss in Sect. [3] where and how the accretion system may form such a beam. 
Finally, in Sect. IHwe will summarize the findings and discuss the strengths and weaknesses 
of the model. 



Electromagnetic model s have recently been discu ssed bv iLevinsoru (120001 ) 



Maraschi k Tavecchiol (120031 ). iKundt fc Krishnal mm . and iKatzl mm . The work de- 
scribed in the following is new in that it starts with simultaneously modeling the X-ray and 
TeV 7-ray data. The data constrain the beam properties tightly, and it becomes possible 
to perform a targeted study of how and where the beam originates. We will describe the 
parallel beams in the AGN frame and we will focus on electron and positrons as emitting 
particles. To make the description more concise, we will frequently use the term electrons for 
both, electrons and positrons. The best studied BL Lac objects Mrk 421, Mrk 501, and lES 
1959-1-650 have all redshifts well below 0.1, and we will omit all redshift and fc-correction 
factors to avoid unnecessary clutter in the equations. 



2. Beam Properties 
2.1. Model Geometries 



We will model the X-ray and TeV 7-ray observations of Mrk 501 taken on April 16th 
1997, the day of the strongest fiare, observed during a spectacular six raonths long fiar- 



ing p e riod. Over the six mont h period several X-ray (jPian et a. 



1997 




Krawczvnski, et al. 




2000[) and 


TeV 


Djannati-Atai et al. 


I999I; 


Ouinn et al. 


1999) 



1998 



Catanese et al. 



1999 



the BeppoSAX (Satellite per Astronomia X) and CAT (Cerenkov Array at Themis) experi- 
ments that observed the April 16 fiare with partial tempo ral overlap. The Mrk 501 data taken 



in 1997 have been modeled by many different grou ps (e.g. lPian et al.lll998l : lKrawczynski et al. 



2OO2I : lMannheinj|l998l : iMiicke &: Protherodl200ll ). We consider two model geometries of the 
emission zones to get an estimate of the model dependencies of the relevant quantities. In 
both cases, we assume a magnetic field configuration with azimuthal symmetry with regards 
to the jet axis. The charged particles move along the magnetic field lines that are nearly 
parallel to the jet axis (the 2;- axis). At the outer edges of the emission zone, the field lines 
make an angle amax to the jet axis. 

In Model 1 the high-energy electrons fill a shape that resembles two back-to-back cones 
with a maximum cone radius of R and a total length of 2 i? along the jet axis (Fig. [1], left 
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side). The specific geometry chosen here reproduces the triangle shaped flares commonly 
observed in X-ray and TeV-7-ray light curves of blazars. Note that the exact shape of the 
region is not very important. The important property of Model 1 is that the diameter divided 
by the speed of light equals approximately the observed flare duration. For simplicity we 
choose an identical height and length of the emitting volume. More realistic shapes may be 
obtained from modeling how an /e~ spark forms in a high electric fleld region. We assume 
that the leptons follow the magnetic fleld lines dissipating little energy until they enter the 
radiation zone at a distance zi from the black hole in which some instability causes the 
particles to move with an isotropic distribution of pitch angles 9 to the magnetic fleld lines 
with 9 < ^^max- We assume that the electrons stop emitting when they exit the radiation 
zone at a distance Z2 from the black hole with Z2 — zi = 2 R. The electrons may stop 
emitting because further instabilities slow them down. The duration Atobs of a flare is given 
by the length of the emission zone: Atobs ~ {^2 ~ zi) / c = 2 R / c. For Mrk 501, typical 
flare durations are Atobs = 12 hrs and thus R ~ 6.5 (Atobs/12 hrs) 10^^ cm. Each electron 
spends a time Atrad = Atobs in the emission region. Here and in the following model, the 
emitting particles travel together with the emitted photons down the jet with the main 
velocity dispersion originating from the pitch angle distribution of the leptons and photons. 

Model 2 describes a spatial electron/positron distribution that might result from a 
synchrotron/pair-creation cascade of extremely high-energy particles accelerated by strong 
electric flelds close to the black hole. We assume that the electrons move along magnetic 
fleld lines that make angles a up to ctmax to the jet axis and are concentrated in a spherical or 
conical shell or "shower front" that travels down the jet axis (Fig. [H right side). We assume 
that the electrons only emit in a radiation zone that extends from a distance di to d2 from 
the black hole. Within the emission zone the electrons spiral with pitch angles 9 < ^max 
around magnetic fleld lines. 

We assume di = ^2/2 and consider only the case where the jet points exactly at the 
observer. High-Lorentz factor electrons emit synchrotron and inverse Compton emission 
only along the direction of their motion. Taking into account the electron pitch angle dis- 
tribution, observers can only see the emission from electrons traveling along fleld lines with 
a < 9ruF,^. Based on the standar d equations from the description of superluminal motion 



(IBlandford. McKee. &: Reed 119771 ) and neglecting (cos(^max)) factors, we get 

K,d2 + T 



Atobs ~ . with 



9r, ^ ^ 



As the electrons travel from di to d2, their different pitch angles result in an increase of the 
thickness of the shell by AT ^ k {d2 — di) . Taking into account that the shell has already 
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a finite thickness at di, we use in simple approximation a constant shell thickness of 

T = (3) 

Combining eqs. ([1]) and ([3j), we get 

max \ / ,\ 

02 ~ ~ 10 — ; — cm (4) 

2 k V12 hrsy V 2° y 

The thickness is T = cAtobs/2 = G-S-IO"*^^ (Atobs/12 hrs) cm. At d2 the shell has the 
following height perpendicular to the jet axis 

cAWsin^^.. ^ f^U^V'cm (5) 

2k V12 hrs J \ 2° J ^ ^ 

For di =6/2/2, each electron spends a time Atrad = {d2—di)/c = Atobs / (4 k) in the emission 
region. In Model 2, triangle shaped pulses could result from a particle density that drops at 
a certain distance from the jet axis. 

Models 1 and 2 represent extreme geometries. The true geometry may lie between these 
two extremes. 



2.2. Simple Analytical Considerations 



In this section, we use some simple analytical estimates to derive insights about the 
beam parameters and their scaling behavior with the model parameters. Further below, we 
will see that the inverse Compton emission is mainly emitted in the Klein-Nishina regime 
where the electrons give most of their energy to the scattered photons. The observation of 
~20 TeV photons thus implies the presence of electrons with Lorentz factors 70 = 4-10'' 
and e nergy Ep = 'jorripC^ = 20 TeV. From the BeppoSAX X-ray observations on April 16, 



1997 flPian et al.lll998l : iMassaro et al.ll2004[ l we infer an 0.1 keV-200 keV energy fiux of 4 = 
4 ■ 10~^ ergs cm~^ sec~^. The observations did not cover the entire fiare, and we assume 
here that the flux averaged over the flare duration was half the peak flux = Ix/2. If the 
source at distance D emits into the solid angle AQ 2 vr (1 — cos V^max 
averaged X-ray luminosity during the flare is 



the time 



Lx = i^Ano^ 



3-10^^ erg sec ^ 



(6) 



The total energy emitted into the X-ray band is 



Wx = Atobs^x 



(7) 



-9- 



Electrons with Lorentz factor 7 moving with pitch angles < ^ < ^max <S 1 emit synchrotron 
radiation at a mean critical frequency 

e^'^ B6 

J ±j "max /oN 

= -1: 8) 

znmeC 

For a Lorentz factor 70, the frequency of synchrotron emission equals Vgo = 2.4-10^^ Hz for 
a magnetic field 

" = Tvft;:;r " ( 2.4 . lo-^ hz ) (rw) (^j « 

where /i is Planck's constant. Averaged over pitch angles, the synchrotron power per electron 
is 



Ps 



The synchrotron cooling time tg = E^/ps is 

In Model 1, the electrons spend a time A^obs in the emission zone and radiate only ~ tg/ 
A^obs ~ 0.05% of their energy before leaving it. Given that the total synchrotron and 
inverse Compton luminosities are approximately equal, we see that the emitting particles 
radiate away only 0.1% of their energy. Model 1 is thus radiatively very inefficient. The 
efficiency does not depend strongly on details of the shape of the emitting region, as long as 
the length of the emission region equals approximately the flare duration. The scaling with 
the model input parameters shows that the efficiency is independent of the electron pitch 
angle distribution. It is lower if electrons with energies exceeding 20 TeV emit the observed 
synchrotron emission. 

In Model 2, the emission zone is more extended. For this model, a rough estimate 
shows that the combined synchrotron and inverse Compton cooling time (ts~^ + ^ic~^) ^ 
equals approximately the time Atobs that the emitting particles stay in the emission 
region. Model 2 is thus radiatively very efficient. 

A total number of electrons 

= , = ^ ^8.5- 10^« (Model 1) (12) 

= 4 K, — ^ 2 • 10^^ (Model 2) (13) 
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produce the flare. The models require a particle luminosity averaged over the duration of 
the flare of 

Le = ^ 6.3 • 10^^ erg sec"^ (Model 1) (14) 



At 



rad 



1.5 ■ 10^^ erg sec"^ (Model 2) (15) 



For Model 2 it is indeed the observed flare duration and not the time At^ad that is relevant for 
computing the power required for sustaining continued flaring activity. These luminosities 
are minimum luminosities to produce the observed X-ray emission. There may be additional 
electrons that produce X-ray emission outside the BeppoSAX energy band. 

The electron beam of Model 2 resembles in some aspects the blobs filled with emitting 
particles of the reference model. This is not too surprising as both models assume that the 
radiation is emitted as synchrotron-Compton emission. The different models result in some 
differences in the phase space distribution of the electrons and photons. Furthermore, the 
beams in Models 1 and 2 could be made entirely of high-energy particles. The reference model 
assumes the presence of a support medium of rather cold particles, and shocks are invoked 
to explain how the energy of the bulk of cold particles is transferred to a few high-energy 
particles. 

It is instructive to compare these jet luminosities to the Eddington luminosity. Based 
on bulge stellar disper sion measurements, the best estimate of the black hole mass i n Mrk 



501 is Mbh ~ 10^ Mq flFalomollKotilainen fc Treves 2002 : iBarth. Ho. fc Sargentll2003h . The 



Schwarzschild radius thus is 



rsch = 2GMbh/c2 = 3-10'' ^^^^ 
and the expression for the Eddington luminosity reads 

LEdd = = 1-25 ■ 10*' erg sec (17) 

o"t IO^Mq 

It is roughly 20 and 10^ higher than the minimum luminosities required by Models 1 and 2, 
respectively. 



2.3. Numerical Results 



While the analytic estimates allow us to understand the scaling of the power- 
requirements, they are not appropriate for describing the inverse Compton component and 
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internal 7-ray absorption processes as these two depend on the details of the synchrotron pho- 
ton energy spectrum. The following numerical estimates use electron energy spectra dNg/d'y 
instead of a mono-energetic electron distribution. We use energy spectra that resemble bro- 
ken power laws over small dynamic ranges with a ratio of the maximum to minimum Lorentz 
factor of ~ 100. The synchrotron power emitted by th e leptons at frequency u pe r frequency 



interval du is computed using the standard equation (IRybicki fc Lightmanlll979l ): 



r^'" , Z"^"'" n ,n dN^ 3e^B sine ^ 
PJu) = ci dj sine de — F(x) 18 

where the first integral runs over the electron Lorentz factors and the second averages over 
the pitch angle distribution. Here and in the following we use the constant ci = (1 — 
cos^max)"^ to normalize the integrals over the pitch angle distribution. The function F{x) 
equals x -^§(0 d^ with the modified Bessel function of 5/2 order, and x = u / Vc- 
The emitted inverse Compton power is approximately given by 

-Pic(i^) = ci / ^7 / dz/s / sine de c [l - cose) a]^^{y)ns{i>s) /\E{y) (19) 

Here we use e'^^^ = ^max in rough approximation. The last term in the integrand is the 
energy that a photon gains in a scattering. The other terms give the scattering rate that 
depends on the angle between the electron and photon velocity vectors and on the density 
of synchrotron photons Ug. The Klein-Nishina cross section is to good approximation 

3 In (2?/) + 0.5 
8 y 

where ctt is the Thomson cross section. The value y is the photon energy in the electron 
rest frame in units of the electron rest mass 

/i z/s 70 (1 -cos^) 

y = 2 21 



Following iDermer &: Schlickeiserl (119931 ) . we use the approximation 



AE{y) = y •y rric for y < land (22) 
= 7 rrie (? for y > 1 (23) 

For Model 1, 11^ can be computed as follows. The total number of photons emitted in a 
certain frequency interval over the duration of the flare is Nslu) = Atobs-Ps(^) / hv. The 
emitting particles occupy a volume of ttR^. Taking into account that the photon density 
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rises from to its final value as the emitting particles move through the emission region, the 
time averaged photon density is 

In the case of Model 2, a similar argument gives after some arithmetic: 



ns(z/) 



2(1 -In 2) AtobsPs(^) 



h V 



(Model 2) 



(25) 



The factor (1 — In 2) arises from averaging the synchrotron photon density over the time the 
emitting particles travel from to d^^ taking into account that the shell height increases 
from if/2 to H and the number of synchrotron photons increases linearly from to its final 
value. Finally, the time averaged fluxes received at Earth can be computed from: 



(Model 1) 



At,ad(Ps(^)+Pic(^)) 

Atobs Af] 



(Model 2) 



(26) 
(27) 



The op tical depth per pa t h leng th for 77 
equations of iGould fc Schrederl (119671 ) 



e+ e processes is computed with the 



Cl 



sin QdQ 



dv a^i^ ns{i') (1 



cos ( 



(2^ 



'0 ^l^thr 

The first integral runs over the pitch angle distribution and the second integrates over the 
target photon frequencies. The threshold frequency for pair creation reads 



and the pair-creation cross section is 

3 (Jt 



— cos( 



(J. 



77 



16 



0') 



2/3(/3'-2) + (3-/3^)/n 



l+(3 
1-/? 



cm 



(29) 



(30) 



with /3 = a/1 — Utxir/y. In the following, we focus only on absorption by synchrotron photons. 
Absorption by external (e.g. disk) photons will briefly be discussed in Sect. 3.2. For Model 
1, the optical depth is approximately Rdri^t/dz and for Model 2 it is ^{d2 — di) dri^t/dz. 

The optical depth for extragalactic pair-creation pr ocesses is computed based on the 
model of the extragalactic background light (EBL) from IKneiske. Mannheim. &: Hartmann 
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( 120021 ) and iKneiske. Bretz. Mannheim. &: HartmannI (120041 ) . This model is almost identical 
to the model P045 of lAharonian et al. ( 20061 ) and is thus consistent with the observed energy 



spectra of two distant BL Lac objects detected with the H.E.S.S. experiment. An equation 
similar to eq. (I5U1) is used with the modifications that the integral over pitch angles goes 
from to TT and the EBL target photon density is used. 

The models together with the BeppoSAX and CAT energy spectra are shown in Fig. 
[21 We use the Beppo SAX data reanalyzed by G . Fossati (private communication, 2006) and 



the CAT data from (IDjannati-Atai et al.lll999l ). In the case of the CAT data, we show the 



systematic (rather than statistical) errors which seem to be the dominant ones. Comparing 
TeV 7-ray energy spectra taken with different Cerenkov telescope experiments (including 
CAT) we obtain a one sigma systematic error on TeV spectral indices of approximately 0.2. 

The difference of the solid and dotted lines in Fig. [2] shows the effect of extragalactic 
7-ray absorption. Although the emitted inverse Compton components peak at and above >5 
TeV, the absorbed energy spectra are rather soft compared to the measured SEDs. Models 
with harder 7-ray energy spectra can be produced with inverse Compton processes deeper in 
the Klein-Nishina regime. We did not implement this possibility here as the resulting models 
do not fully account for the < 1 keV BeppoSAX data, and would thus require additional 
X-ray emission from downstream plasma to reproduce the < 1 keV flux. 

The specific choices of model parameters result in mechanical luminosities of 2.3-10^^ 
erg sec~^ and 2.7- lO^'* erg sec~^ for Model 1 and 2, respectively. For Model 1, we fitted the 
data using Atobs = 22 hrs, 6'max = 1°, «max = 1° as the main input parameters. For Model 
2, we used Atobs = 12 hrs, 6'max = 2°, and Omax = 3°. Compared to the parameters of Model 
1, we had to assume for Model 2 a larger solid angle into which the radiation is emitted and 
a shorter Atobs in order to reproduce the TeV 7-ray flux level. The power requirement for 
Model 2 could be reduced to a value close to the one in eq. (ITBi) with smaller 6'max and Omax 
values and assuming the presence of exter nal seed photons, for example from a n out er slower 



layer of the jet similar as in the model of iGhisellini. Tavecchio. &: Chiabergd (120051 ). 



The main conclusions from the more detailed modeling are that (i) parallel beam syn- 
chrotron self-Compton models are viable without requiring external seed photons to account 
for the observed 7-ray emission, (ii) internal absorption effects are small or negligible, and 
(iii) that the high-energy beams require so much power that power efficient beam formation 
models are strongly preferred. 
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3. Origin of the parallel electron beam 

3.1. Constraints from the Total Energetics 

We concentrate here on electromagnetic models as they predict powerful large-scale 
magnetic and electric fields that might be able to accelerate particles to very high ener- 
gies and produce large-scale ordered motion. Three geometries for accelerating particles in 
electromagnetic models have been discussed in the literature. The models assume that the 
black hole and accretion disk are embedded in a conducting magnetosphere. Free charge 
carriers are created in the magnetosphere through cascades involving curvature radiation or 
inverse Compton scattering and pair creation processes. A quasi-stationary configuration 
is achieved if the charge carriers are distributed such that E ■ B =0. If this condition is 
fulfilled, charged particles move along magnetic field lines without gaining or loosing energy. 
The magnetosphere can act as a series of nearly parallel conductors with zero resistance 
along the magnetic field lines that can sustain a voltage drop far away fro m where i t was 



generated. The voltage drop may either be generated by the disk itself (jLovelacd Il976 



Blandf or dl 1 1 9 76l ) . or by the Kerr black hole spinning in the horizon-thr eading magnetic field 



supported by the disk (IBlandford fc Znajek 



1977 



Thorne et al.lll986h . Katz (2006) argues 



that the electric field generated by a homopolar generator in a rotating magnetized fiuid has 
curls in the observers frame and cannot be shorted out at all locations by a stationary charge 
distribution. It thus seems likely that non-stationary vacuum gaps form, which are capable 
of accelerating particles. A possible geometry of the magnetosphere and the gap region is 
shown in Figure [21 The formation of the gap close to the rotation axis might explain the 
initial collimation of the jet. 

In electromagnetic models, the Poynting fiux launches the jet and anchors the jet to 
the disk or black hole. For simplicity we focus here on disk models; some considerations are 
applicable to the Blandf ord- Znajek model as well. We assume that the coordinate system r, 
0, z is aligned with the jet and the black hole resides at its origin. Directly above the disk, 
the Poynting fiux is 

S=—ExB = —EpXBt (31) 

4 TT 4 TT 

where the subscripts p and t denote the poloidal (r and z) and toroidal (0) components, 
respectively. The second equality follows from the fact that the toroidal electric field vanishes 
for a stationary axisymmetric solution. The condition that the electric field E' in the frame 
co-rotating with the disk material vanishes gives 

= E' = Ep X X Bp (32) 

— * 

with n = (0,0,^2) the angular velocity of the disk material at location r. Thus, the poloidal 
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electric field satisfies ^ 

= -- (n X X Bp (33) 

Combining eqs. (PTi) and the magnitude of the Poynting flux is 

S=^nrBpBt (34) 
47r 



If we now use 



and assume 



O.H^V'^ (35) 



5t(r) = 5t,,, (r/ri)-i (36) 
5p(r) = fip,., (r/n)-^ (37) 

the power transported by the Poynting flux perpendicular to the disk surface is 



^disk 



where we used ri ^ 3-Rsch and r2 ^ ri. Thus, magnetic fields between 10^ and 10^ G 
are needed to form the beam of Model 2. Here and in the following, we use the rather 
high power requirement from the numerical SSC calculations. The reader should keep in 
mind that external Compton models would require ~20 times less power. For Model 1, Bt 
and Bp have to be both roughly ten time s higher. Various estimates o f the magnetic field 



in accretion disks have been discussed by iGhosh fc Abramowiczl (119971 ). Magnetic fields of 
~ 10^ G seem likely from dimensional arguments and numerical simulations. 

The voltage drop across the disk from ri to r2 ^ ri is 



r-2 



1...^ V2iGMBRBp^r, ,^19 ^BH Bp^,, 



V,2= -rQBpdr= ^---^^-v^^ ^ 4.3 . iqi^ ^lEilL y (39) 

c P c2 109 Mq 200 G ^ ' 

If r2 = 2ri (rather than r2 ^ ri), L^^^k and V12 are smaller by a factor of (1 — l/\/2). 
A natural mechanism for explaining the variable nature of the X-ray and TeV gamma- 
ray emission from blazars is screening of the voltage V12 by electron/positron pairs. The 
Blandford-Znajek mechanism results in a qualitatively similar relations between the magnetic 
field, the electric fields, and the emitted power. The main difference is that the EMF is 
generated close to the event horizon rather than at the inner region of the accretion disk. 

Given a poloidal magnetic field with -Bp,ri ~ 200 G at the base of the jet, we can 
estimate the magnetic field at the distance 62 in the simple case that the poloidal magnetic 
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field scales inversely proportional to the radius of the emission zone squared. For Model 2 
we obtain 

where H is the radius of the emission zone perpendicular to the jet axis. This simple estimate 
agrees well with the value inferred from modeling the data (see eq. (|9])). 



3.2. Beam Formation 

The number of electrons required to explain the X-ray and TeV 7-ray emission is so 
large that charge neutrality of the beam is important as otherwise the jet would expand too 
rapidly in the direction perpendicular to the jet axis. The maximum energy to which particles 
are accelerated depend on the relative magnitude of the energy gain and energy loss rates. 
The energy gain rate depends on Vu and the distance over which the voltage drop occurs. 
Energy loss mechanisms are curvature radiation, synchrotron emission and inverse Compton 
emission. With curvature radii on the order of ri, curvature losses are only important at 
very high 3>TeV energies. The magnitudes of the synchrotron and inverse Compton losses 
are highly uncertain as they depend on the magnetic field strength and particle pitch angle 
distribution, and on the intensity of the ambient photon field in the acceleration region, 
respectively. Synchrotron losses are negligible if electrons drift along magnetic field lines. 
Assuming scattering in the Thomson regime, a 20 TeV electrons loses the energy 

AEic = 70^ (1 - COS^Y /acc (Tt Ua (41) 

when traveling a distance /acc through an ambient photon field with energy density Ma. Here, 
i) is the angle between the electron and photon velocity vectors. The energy losses are 
negligible when u^, satisfies 

2 

^ -T, '-l^r - 0-002 erg cm'^ (42) 

70 (1 - COSt')^/accO-T 

for = 30° and /acc = ^i- The ambient photon energy density corresponds to a luminosity 
of 

La = 47r (30 rO^c Ma ~ 5- lO^^erg sec"^ (43) 

if the particle acceleration region is at a distance 30 ri from the photon source, presumably 
the accretion disk. Comparing the minimum beam luminosity required for producing the 
observed X-ray and TeV 7-ray emission (eq. ( ITSl) ) with La, we see that the model requires an 
accretion flow with a radiative efficiency of 3% or less. Higher disk luminosities are possible 
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if the emission frequency is sufficiently high that the high-energy electrons interact only in 
the Klein-Nishina regime. 



Lovelacd (119761 ) assumed that protons may be accelerated all the way to ultra high 
energies and that quasars thus may be accelerators of ultra hig h energy cosmic rays (see 
Boldt fc GhoshI Il999l : iBoldt fc LoewensteinI 120001 : iLevinsonI |2000| for similar recent papers) . 
He stipulated that the flow of high-energy protons may entrain or pick up electrons. As- 
suming that the high-energy electrons and protons would move with identical velocities and 
Lorentz factors, the acceleration of protons would increase the minimum beam luminosity 
by the proton to electron mass ratio. For Model 1, the required luminosity would exceed 
the Eddington luminosity by at least two orders of magnitude. For Model 2, accretion with 
a few times the Eddington rate would be sufficient. With -Bp,ri ~0.17 G, eq. fl5^ predicts 
proton energies of ~ 3 ■ 10^^ eV and just the right electron energy of ~20 TeV. However, a 
prohibitively strong toroidal magnetic field exceeding 10^ G would be required so that the 
Poynting flux can power the massive electron-proton beam. 

The acceleration of electrons or positrons with subsequent entrainment of oppositely 
charged leptons would result in a beam with a much lower power. However, eq. fl39l) predicts 
an adequate voltage drop for a very weak poloidal magnetic field with -Bp^n ~ 2 ■ 10~^ G. 
Even for Model 2, the beam power requires again a toroidal magnetic field exceeding 10^ G. 
Stronger Bp ^ and weaker i?t,ri would be viable if the leptons are accelerated to energies 
exceeding 40 TeV, and then entrain both electrons and positrons, slowing them down to a 
mean energy of ~20 TeV per lepton. 

Acceleration of electrons or positrons with Bt ~ Bp would produce a few particles with 
very high energies. A natural way of transferring the energy from a few high-energy particles 
to many low-e nergy particles are casca d es. Electromagnetic cascades in AGN jets have been 
discu ssed by ( iBurns fc Lovelacd Il982l : iBlandford fc Levinsoru Il995l : iLevinson fc Blandford 
19951). A generic d iscussion of electromagnetic cascades in the >TeV regime has been given 
in (ISvenssonlll987l ). Unfortunately, cascade models need considerable fine tuning to produce 
the electron beams with the right properties. We briefly go through a specific scenario to 
emphasize some of the relevant difficulties. Electrons are accelerated until the energy gains 
in the electric field E12 = V^i2/Lcc along magnetic field lines equal the energy losses. If 
curvature losses dominate, the energy loss rate reads 

2 e'c7^ 
3 



Pc 



p 



(44) 



with p an average curvature radius, which gives a maximum Lorentz factor of 



7n 



3 {V12/L 
2 e 



P' 



1/4 



10 



11 



12 



1019 V 



1/4 



O.lSri 



-1/4 



1/2 



(45) 
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which corresponds to an energy of 50 PeV. The highest-energy electrons emit curvature 
photons at frequency 



An p yiO^^Yj \0.15riJ \ri 

corresponding to a photon energy of ~40 TeV. Overall, each lepton drifting through the 
acceleration region emits 7-rays with a total energy of eVu — 7max'^eC^, while it escapes 
with a relatively small amount of kinetic energy ~ 7max^eC^- A part of the energy in the 
7-ray beam may be converted back to the leptonic sector if the optical depth for pair creation 
is approximately unity. The pair creation cross section (eq. fl30l) ) reaches its maximum value 
of ~ 1/4 (Tt for target photons of frequency z/a ~ 2z/thr- Assuming that the target photons 
have the frequency u ~ 2z/thr, and that the curvature 7-rays are emitted at a distance 30 ri 
from the target photon source, we infer a minimum target photon luminosity of 1.4-10^^ 
erg sec~^ for which the pair creation optical depth for 40 TeV 7-rays escaping to infinity 
equals unity. The model requires fine-tuning as seed photons of frequency ~ 2 i/thr are needed 
to suppress inverse Compton processes of >20 TeV electrons. If the target photon field is 
sufficiently intense, a cascade with several generations of photons and pairs can be initiated. 
The end-product of the cascade then depends on the energy spectrum and spatial gradient 
of the target photons. 

Proton induced s ynchrotron /pair-creation cascades (PIC) in bl azars have been discussed 



by (IMannheiml Il993l : iMiicke &: Protherod I2OOII : iMiicke et al.ll2003l ). These models produce 



7-rays as synchrotron emission of high-energy electrons/positrons with possible contributions 
from other secondary cascade particles. The models require not only high-energy protons, 
but also co-spatially accelerated electrons. The latter emit the radiation field that causes the 
protons to photo-produce mesons, and explain the observed X-ray emission. PIC models had 
originally been proposed in the framework of the shock-acceleration picture. A re-evaluation 
of these models with regards to accelerating the protons and electrons in the strong electric 
fields in the surrounding of a black hole may be a worthwhile enterprise, but is outside of 
the scope of this paper. 



4. Discussion 

This paper discusses the possibility that the X-ray and TeV 7-ray radiation from BL Lac 
type objects is emitted by parallel electron-positron beams that are accelerated by strong 
electric fields close to the accreting central black hole. For the first time we attempt to 
explain both the X-ray and TeV 7-ray emission from these objects with such a model. Fitting 
the model to data from the BL Lac object Mrk 501, we find that the particle acceleration 
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zone and the X-ray and TeV 7-ray emission zone have to be spatially distinct. Otherwise, 
the large magnetic fields required to launch the jets prohibit the simultaneous emission of 
the X-rays and TeV 7-rays as synchrotron and inverse Compton radiation, respectively. 
Previous papers discussing the high-energy emission from blazars dismissed the possibility 
that the high-energy particles producing the observed radiation may be accelerated close to 
the black hole. It was believed that the co-spatially emitted radio to optical emission would 
cause an inverse Compton catastrophe, causing the high-energy particles to lose all their 
energy. For several reasons this argument is not valid for all blazars. First, in contrast to 
the situation in powerful quasars, the accretion disks of BL Lac objects are not radiating 
efficiently. For most BL Lac objects, there are only upper limits on the disk luminosity. 
Furthermore, a common assumption had been that the observed radio to optical radiation 
from BL Lac objects was emitted co-spatially with the X-ray and TeV 7-ray emission. For 
some objects like Mrk 501, this assumption had to be dropped even for the reference model, 
as the upper limits on the size of the emission region from the observed fiux variability time 
scales resulted in synchrotron self-absorption cut-offs in the 100-1000 GHz range. Thus, even 
the standard model requires that at least the radio emission is produced further downstream 
in the jet than the X-ray and TeV 7-ray emission. In BL Lac objects, the observation 
campaigns carried through so far have not produced solid evidence for a correlation between 
the X-ray/TcV gamma-ray fluxes and the infrared-optical fluxes. Thus, also the infrared 
and optical emission may be produced downstream of the X-ray and TeV 7-ray emission, 
or even by independent processes. Two other factors contribute to the suppression of an 
inverse Compton catastrophe. In the parallel beam model the emitting particles and the 
photons travel almost in parallel, greatly reducing their interaction rate; furthermore, inverse 
Compton interactions of >TeV electrons and positrons with all photons with wavelengths 
shorter than a few microns are Klein- Nishina suppressed. 

In the reference model, the jet is launched clcctromagnctically and the jet power is 
subsequently transferred to particles which carry it to large distances from the accretion 
system. Shocks in the jet medium transfer the energy transported by a large number of 
cold particles to high-energy particles that emit the observed radiation. Compared to this 
reference model, the model discussed here avoids the need for transferring the power first 
from Poynting flux to a large number of particles and then to transfer it back to a few high- 
energy particles. Furthermore, the parallel electron-positron beam model can cope with 
some of the difficulties of the reference model: 

1. The model can explain the "odd" energy spectra of emitting electrons/positrons in- 
ferred from synchrotron- Compton fits to BL Lac data with more ease than the Fermi 
acceleration mechanisms. The most striking features are particle distribution func- 
tions almost resembling "delta-functions" or MaxweUian distributions. Such particle 
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distribution might result from the acceleration and/or cascading processes discussed 
above. 

2. The parallel beam model accounts for the non-detection of one of the tell-tale signatures 
of Fermi-type acceleration mechanisms, namely a particularly soft energy spectrum 
during the early rising phases of flares. 

3. A consequence of the nearly parallel flow of the emitting particles and photons in 
the emission region are nearly simultaneous variations of the synchrotron and inverse 
Compton fluxes if fluxes emitted by electrons of similar energies are sampled, and if 
complications arising from cooling of electrons and thus an increase in seed photons 
can be neglected. Synchrotron self-Compton models in which the emitting particles 
move isotropically in the jet frame, pre dict that the synchrotron fluxes should rise 
faster than the inverse Compton fluxes (jCoppi fc Aharonianlll999l ). Although such a 
time lag has long been searched for, it has eluded detection so far. 

4. Sources like Mrk 501 and Mrk 421 show extended flaring periods with many flares. Fur- 
thermore, observations of Mrk 501 in 1997 showed an astonishing stability of the TeV 



7-ray energy spectrum during the entire observation campaign (IKrawczynski. et al. 



2000). The X-ray and TeV 7-ray fluxes followed the same correlation during many 



distinct flares over a time period of several months. In the parallel beam model, regu- 
lar and rather uniform flaring might result from alternating between shorting out and 
evacuating the particle acceleration region. Particle-accelerating recollimation shocks 
at certain typical distances from the central engines provide a viable alternative expla- 
nation. In contrast, mod els in which flares are produced by collisions of plasma blobs 
have difficulties to do so (ITanihata et al.ll2003l ). 



5. In the reference model, the high-energy particles move isotropically. As the parti- 
cle pressure dominates over the magnetic field pressure by rn any orders of magnitude 
( iKino. Takahara. fc Kusunosd 120021 : IKrawczynski et al.ll2002l ). the model has difficul- 
ties to explain why the particles do not flow out of the emission volume with the speed 
of light. In the parallel beam model, the motion of the emitting particles along ordered 
magnetic fleld lines explains the beam coUimation more naturally. 



The parallel beam model has its own challenges. The discussion in the previous section 
indicates that direct particle acceleration followed by entrainment of ambient matter and 
the study of electromagnetic cascades in the TeV-PeV energy regime are areas where more 
detailed modeling is required. Another area for future work is more detailed time resolved 
simulations of the temporal evolution of the beam. 
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In the introduction we mentioned the discrepancy between jet bulk Lorentz factors 
inferred from the reference model to the X-ray and TeV 7-ray data, and from VLBA ob- 
servations. Both, the reference model, and the parallel beam model can solve this problem 
by positing that the emission of the X-rays and TeV 7-rays is the start of a drastic energy 
dissipation of the jet. Model 2 is radiatively very efficient so that a considerable fraction 
of the jet energy goes directly into radiation. If radiatively inefficient models apply, the jet 
may slow down by entraining ambient material. 

In the reference model, the difference between the emission of quasars and BL Lac ob- 
jects is commonly attributed to a difference of the maximum energy of accelerated particles 
owing to more efficient inverse Comptqn coo ling of accelerated particles in the intense radi- 



ation fields of quasars (iGhisellini et al.lll998l ). The parallel beam model has the potential to 
explain not only that the SEDs of quasars and BL Lac objects are different, but also that 
the kpc-scale jets are different. In the reference model, the intense photon fields only affect 
the particles accelerated by the jet far away from the central engine. In the parallel beam 
model, they can affect the process of jet formation itself. Thus, the intense radiation fields 
of quasars may prevent a parallel particle beam from forming altogether, and may allow a 
different jet formation mechanism to produce the powerful quasar jets. 
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Model 1 



Model 2 




Fig. 1. — Two geometries to account for the X-ray and TeV 7-ray emission from BL Lac 
objects. On the left side, the emitting particles fill a volume (shaded grey here) made up 
of two back-to-back cones. The geometry might result from a single larger cone (dotted 
line), where the number of electrons tapers off after some time. The right side shows a 
shell of high-energy particles traveling down the jet. Such a geometry might result from an 
electromagnetic shower developing close to the black hole. In both geometries, the emitting 
electrons and/or positrons drift along the magnetic field lines, and start emitting within the 
regions delimited by Zi and Z2, and by di and ^2, respectively. The true geometry might lie 
between the two extremes shown here. 
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Fig. 2. — Results from modeling the Mrk 501 data from April 16, 1997 with the parallel 
beam model. The left and right panels show fits with Model 1 and Model 2, respectively. 
The solid and dotted lines show the emission with and without accounting for extragalactic 
absorption in 77 e+e" pair creation processes. The April 16 flare was the brightest of 
a large number of flares detected during a 6 month period. We show here the BeppoSAX 
X-ray and CAT '-f-raj fluxes at only 50% of their measured values, to take into account that 
the observations did not cover the entire flare and that the fluxes averaged over the entire 
duration of the flare are likely to be lower than those measured at the peak of the flare. 
In the case of the CAT spectrum the shaded area shows the systematic uncertainties which 
dominate over the statistical uncertainties. The model parameters are given in the text. 
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Black Hole Disk 

Fig. 3. — Snapshot of the magnetosphere and a non-stationary vacuum gap above an ac- 
creting black-hole system. Close to the black hole and the disk, the magnetosphere is filled 
with low-energy e^/e" pairs that move along electric equipotential surfaces (dashed lines). 
Further away, there is a vacuum gap with a voltage drop V12 = V2 — Vi that can acceler- 
ate particles to >TeV energies. See Lovelace, MacAuslan, & Burns (1979) and Lovelace & 
Ruchti (1983) for a slightly different geometry. 



